On the incompatibility of strains and its application to mesoscopic studies of plasticity 
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Structural transitions are invariably affected by lattice distortions. If the body is to remain crack- 
free, the strain field cannot be arbitrary but has to satisfy the Saint- Venant compatibility constraint. 
Equivalently, an incompatibility constraint consistent with the actual dislocation network has to be 
satisfied in media with dislocations. This constraint can be incorporated into strain-based free 
energy functionals to study the influence of dislocations on phase stability. We provide a systematic 
analysis of this constraint in 3D and show how three incompatibility equations accommodate an 
arbitrary dislocation density. This approach allows the internal stress field to be calculated for an 
anisotropic material with spatially inhomogeneous microstructure and distribution of dislocations by 
minimizing the free energy. This is illustrated by calculating the stress field of an edge dislocation and 
comparing it with that of an edge dislocation in an infinite isotropic medium. We outline how this 
procedure can be utilized to study the interaction of plasticity with polarization and magnetization. 

PACS numbers: 05.70.-a, 61.72.Lk, 63.70.+h, 75.85.+t, 81.30.Kf 
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I. INTRODUCTION 

Physical properties of solids are mainly determined 
by crystal defects. Their strain fields couple to the mi- 
crostructure and thus often deteriorate the properties of 
the underlying perfect crystal. The coupling of point 
defects with microstructure can give rise to variation in 
the twin wall widths of ferroelastic martensites . Simi- 
larly, dislocations through their long-range strain fields 
affect the spatial evolution of the microstructure and 
are typically responsible for shifts of the transforma- 
tion temperature . To guarantee phase stability, one 
needs to understand how the phase diagram changes 
as a function of the density of dislocations. This can 
be accomplished efficiently using the continuum theory 
of dislocations in which dislocations are viewed as 
sources of incompatibility of strains and stresses ' . 

The concept of incompatibility of strains was used as a 
basis for the development of field dislocation mechanics 
and continues to serve as a template for the develop- 
ment of mesoscopic models of self-organization of dislo- 
cation networks . It has been readily incorporated 
within the Landau continuum approach to study spin- 
odal decomposition mediated by dislocations as well as 
the effects of dislocations on the stability of phases in 
martensites . These studies have been in two dimensions 
where only one incompatibility relation constrains the 
spatial variation of order parameter fields 16 ' 17 . Analyses 
of the three-dimensional compatibility constraint, which 
guarantees integrability of the strain field in dislocation- 
free media, can be found in . In three dimensions, 
the compatibility constraint is represented by six equa- 
tions whereas only three are required to make the strain 
field integrable . It has not been apparent how this re- 
duction can be justified and this lack of understanding 
has led to conflicting statements about the mutual de- 



pendencies of these equations 22 . Nevertheless, a system- 
atic analysis that examines the structure of the three- 
dimensional incompatibility constraint as applied to a 
finite density of dislocations, is lacking. This has not al- 
lowed for the generalization of two-dimensional models 
of plasticity to three dimensions. The primary objec- 
tive of this paper is to provide a systematic analysis of 
the full set of six incompatibility equations assuming pe- 
riodic boundary conditions and to show how they can 
be reduced to three incompatibility constraints that rep- 
resent an arbitrary dislocation network. This makes it 
possible to develop strain-based free energies 23 ' 24 in 3D 
that describe the spatial evolution of the microstructure 
at the length and time scales that are not amenable to 
atomistic or continuum approaches. 

The elastic free energy due to a distortion from a 
high symmetry parent phase to a lower symmetry prod- 
uct phase is typically written in terms of a finite num- 
ber of order parameters that characterize this change. 
Thus, the transformation from a cubic symmetry to the 
tetragonal, orthorhombic or monoclinic symmetries in- 
volves symmetry- allowed invariants in the free energy, 
the mathematical structure of which is the same for all 
materials that undergo the same type of phase transition. 
Materials undergoing the same crystal symmetry change 
are distinguished by the coefficients of these invariants 
that reflect the temperature and pressure dependencies 
of elastic constants, martensitic strain at the transition 
or curvature of soft phonon modes that drive the tran- 
sition. This is extremely valuable since a single formu- 
lation of the free energy then describes a large class of 
materials. In this paper, we consider the strain-based 
formulation of the Landau- Ginzburg functional for the 
cubic to tetragonal transformation that was developed 
by Barsch and Krumhansl 23 and successfully applied in 
a number of theoretical studies 4 ' 26 ' 27 . 



II. INCOMPATIBILITY CONSTRAINT 

We begin by defining a set of six symmetry-allowed 
combinations of the components of the elastic strain ten- 
sor which can serve as order parameters to characterize 
the symmetry change from the parent (typically cubic) 
phase to an equal or lower symmetry phase. The free 
energy is then written in terms of these order parame- 
ters and its minimization yields the homogeneous ground 
state as well as spatial variations in the order parameter 
fields due to any inhomogeneities, such as gradient terms. 
However, in order to keep the deformed body continuous 
(i.e. without cracks), the components of the total dis- 
tortion tensor f3 t = e t + uj 1 cannot be independent (e* is 
the symmetric elastic-plastic strain tensor and a;* the an- 
tisymmetric tensor of elastic-plastic rotations). Instead, 
they satisfy the equation Vx/3 t = 0, where the total dis- 
tortion is f3 l = (3 J rf3 p . Here, (3 P is the permanent plastic 
distortion of the lattice due to the spatially inhomoge- 
neous distribution of dislocations. It is then customary 
to write V x (3 P = — a, where a is the Nye tensor that 
represents the density of infinitesimal Burgers vectors . 
However, the cohesive forces of the matter cause elastic 
relaxation of the lattice, (3. In order for the total distor- 
tion tensor to be curl- free, it follows that Vx/3 = a. 
This tensor is non-symmetric and its components, o^-, 
correspond to a dislocation whose line direction is par- 
allel to Xi and the Burgers vector parallel to Xj. Hence, 
the diagonal components of a correspond to the three 
screw dislocations with their line directions parallel to 
the three (100) axes, while the six off-diagonal compo- 
nents describe the density of the six variants of edge 
dislocations. Taking the symmetric part of the curl of 
Vx/3 = q yields the so-called incompatibility constraint 



V x V x e = 77 



(i) 



where 77 = sym (V x a) is the incompatibility tensor. 

When developing models of phase transformations 
based on strains, the condition (1) has to be enforced . 
In components, Eq. (1) represents two sets of three equa- 
tions, 
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where z, j, k = {1,2,3} and i ^ j ^ k. The comma means 
a partial derivative and no summation of the repeated 
indices is performed here. In two dimensions, only one 
equation is not satisfied trivially and this represents the 
constraint for the internal strains. The situation is more 
complicated in three dimensions since all six equations 
have to be enforced. However, only three equations are 
needed and the main problem is to find a rigorous way 
to reduce the six equations to three. 

In order to find the minimum number of constraints 
for a spatial variation of the internal elastic strain field 
that complies with a given dislocation density, assuming 
periodic boundary conditions, we begin by writing (2) in 



k-space. In the augmented matrix representation, this 
linear system takes on a block-symmetric form 
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where tilde (~) denotes a Fourier image, and k = 
(fci, &2, ks) is the k-space vector. If no additional con- 
straints on the components fjij are supplied, this system 
of equations would be internally inconsistent. The three 
conditions that make this system consistent are obtained 
from the Bianchi identity V • a = which is only true 
for a permissible dislocation density, i.e. if dislocations 
do not begin or end inside the body. In order to gain 
insight into the structure of (3), let us begin by splitting 
the system into two sets of three equations each, where 
the first set is obtained from (3) by taking rows number 
1-3 and the second set from rows number 4-6. Performing 
the row reduction on each of these subsystems separately 
converts them into the following forms: 



■jyj-gedge _ -edge 
jy/r -screw -screw 

where the matrices of both systems are given by 
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However, the right-hand sides of (4) are different: 
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It is important to emphasize that the right-hand side of 
the first system contains only the diagonal components of 
77. If an edge dislocation is present, only fj edge is nonzero 
whereas f] screw = 0. Similarly, the right-hand side of the 
second system contains only the off-diagonal components 
of r). Hence, for a screw dislocation fj screw is nonzero 
whereas fj edge = 0. Clearly, the first set of equations 
can be used to incorporate edge dislocations, while the 
second set incorporates screw dislocations. However, one 
cannot solve these systems to obtain the strain field in 
terms of the components of the incompatibility tensor r). 
They only represent constraints that the internal strain 
field has to satisfy if a finite density of dislocations is 
present. It is also very important to emphasize that since 
the right-hand sides of these systems are different, the 
constraints imposed on strains by the two systems are in 
general different. Hence, we denote by s edge and s screw 



the total elastic strain fields due to all edge and screw 
dislocations, respectively. Because the two systems are 
linear, the superposition of the two sets yields 
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where e = s edge + £ screw i s the Fourier image of the total 
elastic strain tensor. Eq. (7) represents a linear system of 
three incompatibility constraints whose right-hand side 
contains in general all components of the incompatibility 
tensor. Hence, in three dimensions (7) are three condi- 
tions the internal strain field has to satisfy due to the 
existence of an arbitrary distribution of the dislocation 
density. If these equations are written as constraints, i.e. 
gi[i] = and i = {1, 2, 3}, one can write the free energy 
density for a ferroelastic material in k-space as 



F[i] 




(8) 



where f[i] is the k-space expression of the free energy 
density, A^ the Lagrange multipliers, and Q the volume 
of the k-space. Recall that for two-dimensional problems, 
the sum in (8) contains only one term, namely that corre- 
sponding to fjnj where Xi is the direction of the (straight) 
dislocation line. However, in three dimensions, this sum 
contains three terms that impose the three conditions 
embodied in (7). The development above is completely 
general, applicable to any kind of displacive phase trans- 
formation and an arbitrary (but permissible) dislocation 
density. Below, we use this formulation to outline the de- 
velopment of a self-consistent mesoscopic model for a cu- 
bic to tetragonal transformation in which the microstruc- 
ture is coupled to the density of dislocations. 



III. CUBIC TO TETRAGONAL TRANSITION 

We begin by defining a set of six order parameters 
6i[s] (i = 1,...,6), one of which corresponds to a hy- 
drostatic distortion, two to the change of shape from 
cubic to tetragonal, and three are shear components. 
The concrete expressions for these order parameters can 
be found in 23 . Here, the primary order parameters are 
those that represent the tetragonal distortion of the lat- 
tice, while all others are secondary order parameters. If 
one considers the cubic structure as a reference phase of 
zero energy, then the strain energy density of a distorted 
phase relative to this reference structure can be written as 
f[e] = \cijkiSijSki- Since this is harmonic, its minimum 
corresponds to the reference cubic structure. In order to 
allow for a phase transition to occur, one augments the 
free energy by higher order terms of the primary order 
parameters that respect the cubic symmetry. Finally, the 
length scale is included by adding gradient terms of the 
primary order parameters. The free energy density is 
then written as /[e], where e is a vector of the six order 
parameters. Similarly, we can obtain an equivalent ex- 
pression for the three incompatibility conditions in terms 



of the order parameters. The stationary order parame- 
ter fields corresponding to a fixed dislocation density (i.e. 
incompatibility field 77), are then obtained by variational 
minimization of the free energy equivalent to (8). Since 
the components of the elastic strain tensor are defined in 
terms of the order parameters, the internal elastic strain 
field is straightforward to calculate. The conjugate stress 
field is then obtained at once using the Hooke's law. 

Our objective in the following is to demonstrate how 
the internal strain field is disturbed when a single edge 
dislocation is present. It can be shown that without dis- 
locations, the minimum of the free energy corresponds to 
a twinned microstructure with the normal parallel to the 
(110) direction. Due to the rotational symmetry of the 
parent cubic lattice, we may consider that the normal of 
the martensite plate is parallel to [110]. If we further 
consider a straight edge dislocation whose line direction 
is parallel to £3, all fields are functions of x\ and X2 but 
not of X3. This allows us to reduce the problem to the 
two dimensional plane strain case. 

Let us now insert a [100] edge dislocation into the mid- 
dle of the block, as shown in Fig. la. Here, the Burg- 
ers vector of the dislocation is taken artificially as the 
edge length of one unit cell used to discretize the space 
in our simulation. Due to the periodic boundary con- 
ditions, this represents an array of straight edge dislo- 
cations whose line directions are parallel to xs and the 
Burgers vectors point in the positive x\ direction. Hence, 
the only nonzero component of the Nye tensor is a^i 
and the only nonzero component of the incompatibility 
tensor is 7733. This incompatibility is incorporated into 
the free energy by the procedure explained above. Min- 
imization of the free energy then yields the equilibrium 
internal strain field. This is represented using the or- 
der parameters e^. Since we are interested in the phase 
transformation from the parent cubic phase to the prod- 
uct tetragonal phase, we will be concerned with the spa- 
tial variation of e2, which is proportional to the devia- 
toric strain s\\ — £22- This minimization was first per- 
formed above T c and Fig. la shows the calculated spatial 
variation of the order parameter field e^. If the dislo- 
cation was not present, the minimum of the free energy 
would correspond to the cubic phase (gray color). How- 
ever, if the dislocation is present, the spatially homoge- 
neous microstructure is frustrated by the strain field of 
the dislocation. One can imagine that the dislocation 
was created by sliding into the block an extra half-plane 
(#2, xs) at X2 > which terminates in the middle of the 
block. Hence, the surrounding lattice at x^ > is under 
compression in the x\ direction. Because E\\ < and 
£22 ~ 0, the order parameter 62 becomes negative and 
this corresponds to the blue region in Fig. la. Due to 
the termination of the extra half-plane at X2 = 0, the 
lattice at X2 < is under tension in the x\ direction. 
Since E\\ > and 822 ~ the order parameter 62 is 
positive here and this gives rise to the red part of the 
butterfly shape in Fig. la. If the minimization of the free 
energy was performed below T c and the body was free 
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FIG. 1. Characteristic deviatoric distortion (62) of the lattice caused by an edge dislocation in the middle of the block above 
(a) and below (b) the transformation temperature in a periodic anisotropic medium. As shown in (c), blue (e^, < 0) and red 
(e2 > 0) are the two characteristic tetragonal distortions and gray corresponds to the parent cubic lattice. 



of dislocations, the equilibrium field 62 would represent 
a twinned microstructure in which the two tetragonal 
variants would be separated by twin boundaries of the 
cubic symmetry. As shown in Fig. lb, the presence of 
the dislocation causes a local distortion of this twinned 
microstructure. This results in a misorientation of the 
direction of the tetragonal distortion and an emergence 
of curvature, both in the neighborhood of the disloca- 
tion and also farther away where the tip of the distortion 
meets another twin boundary. In addition to this distor- 
tion, the microstructure accommodates the dislocation 
as a part of the twin boundary. The distortion around 
the dislocation thus represents signatures of the disloca- 
tion in the order parameter field 62 above and below T c . 
For distances large compared to the unit cell in the cal- 
culation, both fields in Figs. la,b approach those for the 
dislocation- free body above and below T c , respectively. 

For completeness, we show in Figs. 2a-c and Figs. 2d- 
f the three components of the calculated internal stress 
field above and below T c , respectively. Here, red repre- 
sents positive stresses, gray the unstressed regions, and 
blue negative stresses. Above T c , the stress field in 
Figs. 2a-c agrees qualitatively with that obtained from 
analytical expressions of the stress field around an edge 
dislocation in an infinite isotropic medium 29 . Below T c , 
the stress field in Figs. 2d-f includes nonlocal interac- 
tions of the the stress field of the dislocation with the 
microstructure. No analytical expression for the stress 
field exists due to the complexity of the problem. As we 
showed already in Fig. lb, an edge dislocation below T c is 
responsible for developing curvature in the system. This 
shows up in the stress field by local bending of the twin 



walls in the vicinity of the dislocation and reversals of 
the state of stress in the direction perpendicular to the 
twin wall. 

Although the calculation above has been performed for 
one edge dislocation, this approach is completely general 
in that the right-hand side of (7) can represent an arbi- 
trary dislocation network. 



IV. GENERALIZATION 

The formalism developed above can be systemati- 
cally extended to study the influence of dislocations on 
the physical properties of multiferroics. To describe 
piezoelectric/piezomagnetic materials in which polariza- 
tion/magnetization is coupled to strain'", one typically 
begins by writing the Gibbs free energy density 

f[Q,e] = /Landau [qf] + 2 ^^Jkl^kl^^^idAA > ( 9 ) 
ijkl i 

where q is a vector of order parameters and /Landau [o\ the 
Landau free energy density consistent with the symme- 
try of the underlying crystal structure. The second term 
in (9) represents the elastic strain energy density and the 
third term incorporates the (in) compatibility constraint 
derived earlier. The spontaneous strain s s is written 
as e s = e — s q , where e is the total elastic strain and 
e ij — QijkiQkQi the strain caused by the change of the 
order parameter. For piezoelectric materials, Q repre- 
sents the electrostrictive tensor and q the polarization 




FIG. 2. Calculated stress field around an edge dislocation in a periodic anisotropic medium (a-c) above and (d-f) below T c . 
Here, (a) and (d) corresponds to an, (b) and (e) to (722, and (c) and (f) to an. The color scale ranges from negative values 
(dark blue), via zero (gray), to positive values (dark red). In all figures the color map spans the same range of magnitudes of 

stresses. 



vector P, while in piezomagnetic materials Q is the mag- 
netostrictive tensor and q the magnetization vector M. 
The first term in (9) then depends only on polarization 
(magnetization), the second term includes a bilinear cou- 
pling of the total elastic strain with polarization (mag- 
netization), and the last term represents a constraint for 
the spatial variation of strain through which dislocations 
affect the order parameter field q. The stationary state 
of the system is thus represented by a generally nontriv- 
ial field of polarization or magnetization, Qi[Qj u^%)-> a ]? 
where a is the Nye tensor of the dislocation density. 



directly obtains the three compatibility equations that 
have to be satisfied at any point in media without dis- 
locations. We have shown how the constraints (7) can 
be systematically incorporated into the strain-based free 
energy functional to obtain a mesoscopic description of 
phase transitions mediated by dislocations. Since the 
lattice constitutes an elastic template that is common 
to all materials, this formalism can be applied to couple 
dislocations not only with microstructure but also with 
other properties of functional materials such as polariza- 
tion and magnetization. 



V. CONCLUSION 
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